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A Quasilinear-Time Algorithm 
for Tiling the Plane Isohedrally 
with a Polyomino 

Stefan Langerman^ Andrew Winslow ^ 


Abstract 

A plane tiling consisting of congruent copies of a shape is isohedral 
provided that for any pair of copies, there exists a symmetry of the 
tiling mapping one copy to the other. We give a 0(n log^ n)-time 
algorithm for deciding if a polyomino with n edges can tile the plane 
isohedrally. This improves on the 0(n^®)-time algorithm of Keating 
and Vince and generalizes recent work by Brlek, Provengal, Fedou, 
and the second author. 


1 Introduction 

The 18th of Hilbert’s 23 famous open problems posed in 1900 1201 concerned 
isohedral tilings of polyliedra where every pair of copies in the tiling has 
a symmetry of the tiling that maps one copy to the other (see Figure [^. 
Hilbert asked for an example of an anisohedral polyhedron that admits a 
tiling, but no isohedral tilings. 

Reinhardt [31] was the hrst to give an example of an anisohedral polyhe¬ 
dron. Along with this example, Reinhardt also stated that a proof that no 
anisohedral polygons exist was forthcoming, a claim thought to be supported 
by Hilbert [T5|. In fact, Reinhardt (and Hilbert?) were mistaken: no such 
proof is possible and Heesch provided the hrst counterexample in 1935 [TH| 
(see Figure [^. 
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Figure 1: Isohedral (left) and anisohedral (right) tilings of a polyomino. 
There is no symmetry of the right tiling mapping one colored tile to the 
other. 

In the 1963, Heesch and Kienzle [19] provided the hrst complete classi- 
hcation of isohedral tilings. This classihcation was given as nine boundary 
criteria: conditions on a polygon’s boundary that, if satished, imply an iso¬ 
hedral tiling and together form a necessary condition for isohedral polygons. 
Each boundary criterion describes a factorization of the boundary into a 
specihc number of intervals with given properties, e.g., an interval is rota- 
tionally symmetric or two intervals are translations of each other. Special 
cases of this classihcation have been rediscovered since, including the crite¬ 
rion of Beauquier and Nivat | 3 | and Conway’s criterion, attributed to John 
H. Conway by Gardner (H] |33] (see Figure [^. 

While a complete classihcation of isohedral tilings exists, many problems 
in tiling classihcation and algorithmics remain open. For instance, com¬ 
plete classihcations of pentagons that tile the plane were claimed as early 
as 1968 [ 23 ], but additional pentagons have been discovered as recently as 
2015 [27]. The existence of an algorithm for deciding if a polyomino tiles the 
plane is a longstanding open problem [laiis], as is the existence of a polygon 
that tiles only without symmetry [35] . 

One of the most successful lines of work in tiling algorithmics was initi¬ 
ated by Wijshoh and van Leeuwen [56], who considered tiling the plane using 
translated copies of a polyomino (isohedrally or otherwise). They proved that 
deciding whether a polyomino admits such a tiling is possible in polynomial 
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Figure 2: The anisohedral polygon of Heesch [TH] and an anisohedral poly- 
omino of Rhoads |32]. There is no symmetry of either tiling mapping one 
colored tile to the other. 

time. Their algorithm was subsequently improved by Beauquier and Ni- 
vat [3], who gave a simple boundary criterion for polyominoes that admit 
such a tiling. Subsequent application of more advanced algorithmic tech¬ 
niques led to a sequence of improved algorithms by Gambini and Vuillon m, 
Provengal [30], Brlek, Provengal, and Fedou [5], and the second author 1371. 
who gave an optimal 0(?7,)-time algorithm, where n is the number of edges 
on the polyomino’s boundary. 

The boundary criterion of Beauquier and Nivat matches one of the cri¬ 
teria of Heesch and Kienzle, implying that this problem is a special case of 
deciding if a polyomino is isohedral. The general problem of isohedrality 
was proved decidable in 1999 by Keating and Vince 121 , who gave a matrix- 
based algorithm running in (9(n^®) time. Their algorithm does not make 
use of boundary criteria, which we note yields a straightforward 0(n®)-time 
algorithm. 

Here we give a 0(?7, log^ n)-time algorithm for deciding if a polyomino 
is isohedral. The algorithm uses the original boundary characterization of 
Heesch and Kienzle [19] to decompose the problem into seven subproblems, 
each of recognizing whether a polyomino’s boundary admits a factorization 
with a specihc form. Structural and algorithmic results on a variety of word 
problems are used, extending the approach of 1371 to factorizations of six 
additional forms. The algorithm also hnds a witness tiling and is easily 
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extended to other classes of lattice shapes, e.g. polyhexes and polyiamonds. 


2 Definitions 

Although the main result of the paper concerns highly geometric tilings, 
the proof is entirely described using words, also called strings. We use the 
term “word” for consistency with terminology in previous work on tilings of 
polyominoes. 

Polyomino and Tiling. A polyomino is a simply connected polygon 
whose edges are unit length and parallel to one of two perpendicular lines. 
Let = {Ti,T 2 ,...} be an inhnite set of hnite simply connected closed 
sets of Provided the elements of ^ have pairwise disjoint interiors and 
cover the Euclidean plane, then is a tiling and the elements of are 



B, C, D, E palindromes 


Figure 3: Two of seven boundary criteria characterizations of isohedral 
tilings. These criteria were given by Beauquier and Nivat |3] (top) and John 
H. Conway [11] (bottom). Precise dehnitions are given in Section]^ 
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called tiles. Provided every Tj G ^ is congruent to a common shape T, then 

is monohedral and T is the prototile of In this case, T is said to 
have a tiling. A monohedral tiling is also isohedral provided, for every pair 
of elements Tj,Tj G there exists a symmetry of ^ that maps Tj to Tj. 
Otherwise the tiling is anisohedral. 

Letter. A letter is a symbol x G S = {u, d, 1 , r} representing the direc¬ 
tions up, down, left and right. The Q°-rotation of a letter x, written te{x), 
is dehned as the letter obtained by rotating x counterclockwise by 0°, e.g., 
O 7 o(u) = r. A special case of ©“-rotations is the eomplement of a letter, 
written x and dehned as x = ti8o(x). The Q°-reflection of a letter x, written 
/©(x), is dehned as the letter obtained by rehecting x across a line with angle 
0° G {— 45°, 0°, 45°, 90°} counterclockwise from the x-axis, e.g., /45(u) = r. 

Word and Boundary Word. A word is a sequence of letters and the 
length of a word W, denoted \W\, is the number of letters in W. For an 
integer i G {1, 2,..., |hF|}, W[i] refers to the Ah letter of W and hF[—z] 
refers to the Ah from the last letter of W. The notation IT* denotes the 
repetition of a word i times. In this work two kinds of words are used: non¬ 
circular words and circular words (dehning the boundaries of polyominoes). 
A word is non-circular if it has a hrst letter, and circular otherwise. For a 
circular word IF, an arbitrary but hxed assignment of the letter IF[1] may 
be used, resulting in a non-circular shift of IF. The boundary word of a 
polyomino P, denoted B{P), is the circular word of letters corresponding 
to the sequence of directions traveled along cell edges during a clockwise 
traversal of the polyomino’s boundary (see Figure]^. 


Figure 4: Polyominoes with (circular) boundary words ur^d^lu^l^, u^r^dl^dl, 
and u^(rd)^ldl (from left to right). 

Rotation, complement, and reflection. The rotation (or comple¬ 
ment, reflection) of a word IF, written ^©(IF) (or IF, /©(IF)), is the word 
obtained by replacing each letter in IF with its rotation (or complement, re¬ 
flection). The reverse of a word IF, written IF, are the letters of IF in reverse 
order. The backtrack of a word IF is denoted IF and defined as IF = IF. 
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Factor. A factor of W is a contiguous sequence X of letters in W, written 
X :< W. For integers I < i,j < \W\ with i < j, W[i..j] denotes the factor 
of W from W[i\ to fF[j], inclusive. A factor X starts or ends at W[i\ if 
W[i] is the hrst or last letter of X, respectively. Two factors X,Y may 
refer the same letters of W or merely have the same letters in common. In 
the former case, X and Y are equal, written X = Y, while in the latter, X 
and Y are congruent, written X = Y. For instance, if IF = uuulruuu then 
IF[1..3] = IF[ 6 .. 8 ]. A factorization of IF is a partition of IF into consecutive 
factors Fi through Fk, written IF = F 1 F 2 ... Fk. 

Prefix, sufRx, affix, middle, and center. A factor X ^ IF is a prefix 
if X starts at IF[1], written X ^pre IF. Similarly, X ^ IF is a suffix if X 
ends at IF[—1], written X IF. A factor X ^ IF that is either a prefix 
or suffix is an affix, written X ^aff IF. A factor X ^ IF that is not an affix 
is a middle, written X ^mid IF. The factor X ^ IF such that IF = UXV, 
\U\ = |F|, and |X| G {1,2} is the center ofW. Similar dehnitions for words 
are defined equivalently, e.g., a word is a prefix of another word provided it 
is congruent to a prefix factor of that word. 

Period, composite, and primitive. A word X is a period of W 
provided IF is congruent to a prefix of X^ for some k > 0 (introduced 
by [21]). Alternatively, X is a prefix of IF and W[i\ = IF[z + |X|] for all 
l<i<|IF| — |X|. A word X is composite provided there exists a subword 
Y such that X = Y^ for some k >2, and otherwise is primitive. 

0-drome, square, and mirror. A word X is a Q-drome provided X = 
Fte+i 8 o(F), e.g., a palindrome is a 180-drome. Such a factor is admissible 
provided IF = XU with f/[—1] 7 ^ te+i 8 o(t/)[l]|^ 

A word X is a square provided X = Y^ for some Y, and a reflect square 
provided X = F/©(F) for some 0. Such a factor is admissible provided 
IF = XF with F[-l] 7 ^ fe{U[-l]) and U[l] ^ feiY)[l]. 

A factor X ^ IF is a mirror, written X IF, provided IF = XUYV 
with Y = X and \U\ = \V\. Such a factor is admissible provided [/[I] 7 ^ 

cTFi], v[i] # vpi] 

Gapped mirror and (reflect) square. A pair of disjoint factors 
X,Y ^ IF is a gapped mirror provided X = Y. Such a pair X, Y is ad¬ 
missible provided IF = XUYV with U[l] ^ U[-l], F[l] ^ F[-l]. 

^In this work, several types of factors and pairs of factors have restricted “admissible” 
versions. Intuitively, admissible versions are maximal in a natural sense for each type. For 
instance, a 0-drome is admissible if it is the longest 0-drome with its center. 
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A pair of disjoint factors X,Y ^ ly is a gapped (reflect) square provided 
X = Y {X = /©(y)). Such a pair is admissible provided W = XUYV with 
(7(1] ^ rll] (C/|l| # /e(V'll])) and C/|-l] ^ V'l-l] ([/[-I] # /e(V'|-l])). 

3 Proof Overview 

The remainder of the paper is dedicated to proving the following theorem: 

Theorem 3.1. Let P be a polyomino with \B{P)\ = n. It can be decided in 
O(nlog^n) time if P has an isohedral tiling (of the plane). 

Here we survey some of the ideas involved in the proof. The proof starts 
with a list of the boundary word factorization forms that together charac¬ 
terize the polyominoes capable of tiling isohedrally. These are found in the 
bordered subregiorj^ of Table 10 of |19j^ excluding the two types of isohedral 
tilings that use 60° and 120° rotations of the shape. The factorizations can 
be cross-verihed using the incidence and adjacency symbols of a more de¬ 
tailed classihcation of isohedral tiling types of Griinbaum and Shephard di, 
and correspond to the isohedral types IH 1, 4, 28, 2, 3, 5, and 6 in this 
classihcation. The factorization forms are: 

• Translation: ABC ABC. 

• Half-turn: W = ABCADE with B, C, D, E palindromes. 

• Quarter-turn: W = ABC with A a palindrome and B, C 90-dromes. 

• Type-1 rehection: W = ABfQ{B)ACf^{C) for some 0, <h. 

• Type-2 rehection: W = ABCAfQ{C)fQ{B). 

• Type-1 half-turn-rehection: W = ABCADf^^D) with B, C palin¬ 
dromes. 

• Type-2 half-turn-rehection: W = ABCDfsi^B) frs>{D) with A, C palin¬ 
dromes and 0° — 4)° = ±90°. 

translation of the caption of Table 10: “The strong border contains 9 major types, 
from which the others can be thought of as emerging by shrinking lines or line pairs.” 
^Reproduced on page 326 of [34] . 
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Winslow [ 37 ] gave a (9(n)-time algorithm for deciding if a bonndary 
word of length n has a translation factorization. Sections throngh give 
0{n log^ n)-time or faster algorithms for deciding if a bonndary word of length 
n has a factorization of each remaining form. Brlek, Koskas, and Provengal |1| 
provide a 0(?7,)-time algorithm for deciding if a given circnlar word is the 
bonndary word of a polyomino, and we assnme for the remainder of the pa¬ 
per that the inpnt is gnaranteed to be the bonndary word of a polyomino 
and thus simple. This assumption of simplicity is used to prove that factors 
and pairs of factors in a factorization are admissible: maximal in a natural 
sense for each factor (pair) type. E.g., for half-turn factorizations: 


Lemma 


4.14 


Let P be a polyomino and B{P) = ABCADE with B, C, D, E 
palindromes. Then the gapped mirror pair A, A and palindromes P, C, D, 
E are admissible. 


For various factors and factor pairs, this implies there are 0{n) or 0{n log n) 
candidate factors for these elements of a factorization, and they can be com¬ 
puted in similar time. Note that this alone is not sufficient to solve these 
problems in quasi-linear time. Without additional structural results, at¬ 
tempting to combine a quasi-linear number of factors into a factorization 
with one of the 6 forms is at least as hard as searching for cycles of hxed 
length between 3 and 6 in a graph with n vertices and 0{n) edges, only 
known to admit a 0(?7,^'®^)-time algorithm pQ. Additional structural results 
must be used, such as the following for half-turn factorizations: 


Lemma 4.4 The prefix palindrome factorization of a word W has the form 


W = with: 


• Every Xi primitive. 

• For all 3 < i < m, ^ 1-^*1 thus m = 0(log \W\). 


Such results allow more efficient “batch processing” of factors to achieve 
quasi-linear running time. It can also be seen by cursory examination that 
each algorithm returns affirmatively only once witness factorization is found. 
Witness factorizations dehne the set of boundary intervals shared by pairs of 
neighboring tiles in an isohedral tiling, thus the algorithm can also return a 
witness isohedral tiling if desired. 




4 Half-Turn Factorizations 


Definition A half-turn factorization of a boundary word W has the form 
W = ABCADE with B, C, D, E palindromes. 

4.1 Prefix palindrome factorizations 

Definition Let W he a word. A factorization W = E 1 E 2 ... En+i is a prefix 
palindrome factorization of W provided that the set of prefix palindromes of 
W is {E 1 E 2 ... Ei : I < i < n}. 

Lemma 4.1. Let W = PX with P, W palindromes and 0 < |P| < \W\. 
Then W has a period of length |X|. Furthermore, if X is composite, then W 
has a prefix palindrome longer than P. 

Proof. Since W = PX anc^P, W are palindromes, W = W = PX = XP = 
XP. So P is a prefix of XP and so X is a period of P and oi W = XP. 
Since X is a period of W and |X| < |1X|, there exist words Y, Z such that 
X = ZY,W = {ZYyZ, and P = {ZYf-^Z for some p > 1. ho X = YZ 
and since W is a palindrome, Y and Z are palindromes. 

If X is composite, then X = for some k >2. So there exist words Gi, 
G 2 such that G = G 1 G 2 , Y = (GiG' 2 )*Gi, Z = G 2 (G'iG 2 )^~*~^- Since Y and 
Z are palindromes, Gi and G 2 are palindromes. 

Now we construct a prefix palindrome of W, called Q, that is longer 
than P. Withont loss of generality, assume IG 2 I > 0. Then there are two 
possibilities for the valnes of IG 2 I and jZj: 

1. IG 2 I < \Z\ and we let Q = X^G 2 . 

2. IG 2 I = jZj and we let Q = PG 2 . 

In the first case, Q = X^G 2 = (GiG 2 )^G '2 = (G' 2 G'i)^^G' 2 , so Q is a 
palindrome. Also, W = {ZYyZ and G 2 is a prefix of Z, so Q is a prefix of 
W and |P| = |IT| - |X| = |IT| - \ZY\ = p\ZY\- |X| < \Q\ = p\ZY\ + IG 2 I < 
p\ZY\ + \Z\ = \W\. 

In the second case, Y = G 2 and Z = SoQ = PG 2 = (ZYy~^ZG 2 = 

and Q is a palindrome. Also, W = PX and G 2 is a 

prefix of X, so Q is a prefix of W and |P| < IQI = |P| + IG 2 I < |P| + |X| = 

|IT|. □ 
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The following is a well-known result; see Chapter 2 of Crochemore and 
Rytter [ 6 ]. 

Lemma 4.2 (Fine and Wilf’s theorem [S]). Let W be a word with periods of 
length p and q. If p + q < \W\, then W also has a period of length gcd(p, q). 

Lemma 4.3. Let Pi,P2, ..., Pm be the set of prefix palindromes of a word 
with 0 < |Pi| < IP 2 I < ••• < \Pm\- Then for any 1 < i < m — 2, either 

l-Rt+ll ~ \Pi\ = |-Pi-|-2| ~ |-Pi+l| or |Pj| -|- |Pii+l| < |Pi!+2|- 


Proof Let Pj+i = P^X^ and Pi +2 = PiXiXi+i. 

Since there are no prehx palindromes of length between |Pj| and l-Pi+ll 
or |Pi+i| and |Pi+ 2 |, Lemma 4.1 implies Xi and Xj+i are primitive and Pj+i 
and Pi +2 have periods of length \Xi\ and |Xi+i|, respectively. Since Pj+i is a 
prehx of Pi+ 2 , it also has a period of length |Xi+i|. 

The lemma permits |Xi| = |Xi+i|, so assume |Xi| 7 ^ |Xi+i|. If |Xi+i| < 

I Pi I, then Pi+i has periods of length \Xi \ and |Xi+i| with |Xi+i|-|-|Xi| < |Pi+i|. 
Then by Lemma 4.2, Pj+i has a period of length gcd(|Xi|, |Xi+i|). This length 
must be at least |Xi| and |Xi+i|, otherwise Xi or Xi+i is not primitive. So 
|Xi| = |Xi+i|. Otherwise, |Xi+i| > |Pi| and so \Pi\ + \Pi+i\ < |Xi+i|-h|Pi+i| = 

IP 


i+2\ 


□ 


The next lemma is a strengthening of similar prior results by Apostolico, 
Breslauer, and Galil [2], I et ah and Matsubara et al. [28] . 

Lemma 4.4 (Prehx Palindrome Factorization Lemma). The prefix palin¬ 
drome factorization of a word W has the form W = X^^X^'^... X^Q with: 


Every Xj primitive. 


• For all 3 < i < m, X]j=i ^ 1^*1 ond thus m = 0(log |IF|). 

Proof. We give a constructive proof. Let Pi, P 2 ,..., P„ be the set of prehx 
palindromes of W with |Pi| < IP 2 I < ■ ■ ■ < |P„|. 

Let Wi be the word such that Pi+i = PilFi and let Q be the word such 
that W = PnQ. So W has a prehx palindrome factorization W 1 W 2 ... WnQ. 


By Lemma 4.1, every IFi is primitive. Moreover, by Lemma 4.3 either |hFi| = 
|IFi+i| or |Pi| < |IFi+i| for every 1 < z < n — 2 . 

Suppose I IFi I = By Lemma 4.3, Pi and Pi+i have a common 

period and thus IFi = IFi+i. More generally, if |IFi| = |IFi+i| = ■■■ = 
|IF,+e|, then W = IFi+i = ••• = IFi+,. If |Pi| < |IFi+i|, then |IF,+i| > 
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\Pi\ = \^j\- factorization • • • ^nQ can be rewritten as 

X{^X2^ ... X^Q with the property that \Xi\ > So for all f > 4, 

2 |Xj_ 3 | < \Xi\ and thus m = 0(log |hh|). □ 

Such a factorization can be stored using 0(log |hh|) space by simply stor¬ 
ing \Xi\ and Ti for each i. Additional observations can be used prove that |fh| 
prehx palindrome factorizations of the suffixes of a word W can be computed 
in optimal time: 

Lemma 4.5. The prefix palindrome factorizations of all shifts of a circular 
word W can he computed in 0(|fh| log \ W\) total time. 

Proof. Lemma 9 of [21] states that the prehx palindrome factorization of 
a non-circular word xY can be computed in 0(log |y|) time given the fac¬ 
torization of Y. Thus the factorizations of non-circular word WW can be 
enumerated in 0(|hL| log \W\) time, beginning with Y = 1T[—1]. Every shift 
of word hh is a subword of the non-circular word WW, and the computed 
factorizations can be trimmed in 0(log |hL|)-time per factorization to be the 
factorizations of shifts of W. □ 

Identical results, including a suffix palindrome factorization lemma, clearly 
hold for suffix palindromes as well. 


4.2 Algorithm 

The main idea is to iterate over all pairs of adjacent letters and guess the 
form of the palindromes D and E in both directions from that location. 
Specihcally, guess what repeated factor XT terminates them in their prehx 
and suffix palindrome factorizations. Then try to complete the factorization 


using Lemma [4.13[ which decides if it is possible to rewrite a given portion of 
the boundary as LfABCAR^ with 5, C palindromes and b, c in some range. 
The results leading up to Lemma 4.13] provide the necessary structure to 
achieve this goal. In particular. Lemma |4.11| shows how to decompose a 
word into two palindromes, and Lemmas |4. 6| through |4.8| yield fast detection 
of a factorization of the form BCR^ with 5, C palindromes. 


Lemma 4.6. Let W he a word with subwords L, R such that W = LR^ and 
R Z^suff L. Let Pi, P 2 be palindromes such that W = PiP 2 R^ with \L\ < |Pi|. 
Then there exists a palindrome P 2 and integer k' such that W = PiP^R^ 
with |Pj| < |P|. 
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Proof. Since |L| < |Pi|, P2 is a a suffix of P* for some minimal i. If i < 1 , 
then either IP2I = |P| (and IP2I = 0 , P = fc + 1 ) or IP2I < |P| and the claim 
is satished. If z > 2 , then there exist words Y, Z with |y| >0 such that 
R = YZ and P 2 = Z{YZy-\ 

Let = Z and k' = i - I + k. So hh = PiP2P^ = PiP^(FP^)*-^P^ = 
PiP^P^'. Since \YZ\ = \YP^\ = |P| and |F| > 0, it follows that |P^| < \R\. 
Since P2 = Z{YZy~'^YZ, Z = P2 is a palindrome. □ 

Lemma 4.7 (Lemma C 4 of [ 9 ]). If a word X1X2 = Y1Y2 = Z1Z2 with X2, 
Yi, Y2, Zi palindromes. Then Xi and Z2 are palindromes. 

Lemma 4 . 8 . Let R be a primitive word and let W = LK^. There is a set of 
integers H with \H\ = 0 (log |IL|) such that W = PiP 2 P^ if and only if it does 
so with \LR^\ < I Pi I < for some h E H. Moreover, given |P| and 

the prefix palindrome factorization ofW, H can be computed in 0 (log \W\) 
time. 


Proof. We may assume IP 2 I < |P| by Lemma 4.6 Consider the prehx palin¬ 
drome factorization of W as described in Lemma 4.4[ Any solution Pi ends 
with one of the repeating subwords Xi of the factorization. There are three 
cases: \Xi\ < |P|, \Xi\ > |P|, and \Xi\ = |P|. 

Case 1: |Xj| < \R\. We claim that if \Xi\ < |P|, then Xy overlaps P*" 
in at most two repetitions of P (and there are at most three values of h). 
Assume, for the sake of contradiction, that P^ is a subword of \Xy\. Then 
P^ is a word of length at least |Xi| -|- |P| with periods of length |X, 


So by Lemma 4.2, P has a period of length gcd(|Xi|, |P|) < \Xi 
contradiction. 


and |P|. 
< |P|, a 


Xil > |P|. We claim that if |Xj| > |P|, then Xi cannot 

Assume, for the sake 


Case 2: 

repeat in P'’ (and there are at most two values of h) 
of contradiction, that Xf is a subword of . Then by Lemma 4.2 


has a period of length gcd(|Xj|, |P|) < |P| < |Xj|. So by Lemma 4.1, the 


factorization given was not a prehx palindrome factorization, a contradiction. 
Case 3: |Xi| = |P|. We claim that if |Xj| = |P|, then h = 0 suffices. 


Suppose |LP^| < 1^11 < ILP^+^I for some h > 1. By Lemma |4.l[ Pi has a 


period of length \Xi\ = |P|. Let Y, Z be words such that YZ vs, a period 
of Pi and |P| — \Z\ = IP2I. So Pi = (YZyY for some p > 1 and Y, Z are 
palindromes. 

Since LR^+^ = {YZ)pYP 2 and \YZYP2\ = 2\R\, YZ = YP2 = R. So 
LR = (YZy~^YP2 = P(P 2 , where P[ = (YZy~^Y and thus is a palindrome. 
So there exists a P[ with |LP*’| < |P(| < |LP^|. 
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Computing H. The value of h in case 3 is always 0. For case 1, use 
the values of h such that Xh contains the last letter of LR^ or For 

case 2, use the values of h such that the prehx palindrome ending at the 
unique repetition of Xi has length between \LR^\ and \LR^^^\. □ 


Lemma 4.9. Let R he a primitive word and let W = LRL. Assume that |-R| 
and the prefix and suffix palindrome factorizations of W are given. Then it 
can be decided in 0(log|hF|) time if W = PiP 2 R^ with Pi, P 2 palindromes 
and \L\ < |Pi|. 


Proof. We may assume IP 2 I < |-R| by Lemma 4.6 First, use Lemma 4.8 to 


compute a 0(log |fF|)-sized candidate set of integers H such that a solution 
exists if and only if = P 1 P 2 with |LP^| < |Pi| < for some 


h & H. By Lemma 4.7, it suffices to check for such solutions with at least 


one of the following types of palindromes: 


• The longest prehx palindrome of with length at least |LP^|. 

• The longest suffix palindrome of with length less than |P|. 

Compute the longest prehx palindromes of LR^~^^ for all values of h in 
0(log|hF|) total time using a two-hnger scan of (1) the prehx palindrome 
factorization of W and (2) the values of h. Use a second two-hnger scan of 
(1) these prehx palindromes and (2) the suffix palindrome factorization of 
the last \R\ letters of the suffix palindrome factorization of LR^ to search for 
a solution Pi, P 2 . 

The longest suffix palindrome of (with length less than R) is in¬ 

variant for h and can be computed in 0(log|hF|) time, using the last |P| 
letters of the suffix palindrome factorization of LR^. Call the length of this 
palindrome A. Use a scan of the prehx palindrome factorization to determine 
if a prehx palindrome of W has length |LP^| — A for some value of h. □ 


Next, we develop a second result that is combined with the previous 
lemma to obtain Lemma 14.131 


Lemma 4.10. Let L and R be words such that L 2 ^pre R. Let Ai he the 
longest common prefix of ffi~'^R and a word U. Let k = I — [lAol/lLl]. Then: 

• For all i with 0 < i < k, Ai = Aq and |^j| < \L^~^\. 

• For all i with k + 2 < i < I, \Ai\ — = |v4fc_|_2| — 
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Proof. Let k be the maximum k such that |Afc| < We show that this 

value of k has the desired properties, including that k = l — [lAol/l-Ll] . 

Property 1. Let 0 < z < fc. Since At is a prehx of V~^R and \Ai\ < 
Ai is the longest common prehx of V' and U. This is true for all 
choices of Ai, and thus all Ai are equal. 

Property 2. By dehnition, |v4fc_|_i| > and so is a prehx 

of U. Since L 2 <pre R, the length of Ak+ 2 ) the longest common prehx of 
]J-{k+ 2 )j^ and U, must be less than So |Afc+ 2 | < and L 

is a period of Ak+ 2 - 

Let Ri, i ?2 be words such that R = R 1 R 2 and Ak +2 = Then 

|/?i| < \L\ and so Ri is the longest common prehx of R and L. 

So for all z > k + 2, the longest common prehx of L^~''R and is 

R~^Ri. Moreover, since and is a prehx of U, 

the longest common prehx of V~'’R and U is also V~''Ri. 

Property 3. Finally, we prove that k = I — [lAol/lLl] . Since l^ol = 
\Ai\ < for all 0 < z < /c, it follows that l^ol < \V‘~^\ = \L\{l — k). Then 
by algebra, k < I — |74o|/|L|. 

Since Ak+i > jt must be that U has a prehx So 

|y4o| = \Ak\ > = |L|(/ —fc —1). Then by algebra k + 1 > l — \Ao\/\L\. 

So k < I — I^dol/ILI < k + 1. □ 


Lemma 4.11. Let W be a word and li, I 2 integers. Assume the prefix and 
suffix palindrome factorizations of W are given. Then it can he decided in 
0(log|hF|) time if there exist palindromes Pi, P 2 such that W = P 1 P 2 with 

l-Pil ^ h, I-P 2 I ^ 12- 


Proof. By Lemma 4T, such a pair of palindromes exist if and only if there 
exists such a pair such that either Pi is the longest prehx palindrome of 
W with I Pi I < I IF I — /2 or P 2 is the longest suffix palindrome of W with 
IP 2 I < |hF| — li. Scan each factorization in 0(log|lF|) time to hnd these 
specihc palindromes, and then scan the opposite factorizations for a second 
palindrome to complete W. □ 


The following result comes from a trivial modihcation of Theorem 9.1.1 
of [T7j to allow for circular words, namely giving the concatenation of two 
copies of a corresponding non-circular word as input, and returning 00 if the 
output has length more than lcm(|X|, |F|). 
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Lemma 4.12 (Theorem 9.1.1 of [IT]). Two circular words X, Y can he 
preprocessed in 0(|X| + |F|) time to support the following queries in 0(1)- 
time: what is the longest common factor of X and Y starting at X[i] and 

YW 

Lemma 4.13. Let W be a circular word. Let W = VZ and W = YR!' such 
that L Z, R :^sus Y, and L, R are primitive. Assume that the prefix 
and suffix palindrome factorizations of every shift of W are given. It can be 
decided in O{{I+ r) log \ W\) time if there exist positive integers b, c such that 
W = Lf APiP2AR^ with A, A admissible and Pi, P2 palindromes. 


Proof. The approach is to iteratively search for a solution for each value 
of b, carrying out the same algorithm on each value. Then performing an 
identical, symmetric iteration through the values of c. First assume fe is a 
hxed value and c is not. Let Ai be the longest word such that L’^Ai and AiR^ 
are a prefix and suffix of W, respectively. Lemma 4.10 implies that there 
exists an integer k such that: 

• For all i with 0 < i < k, is hxed and = l^ol -|- \R\i. 


• For all i with k + 2 < i < I, |v4ji?*| is hxed and |L^y4j| = \L^\ + |Afc +2 
\R\(i-{k + 2}). 


• k can be computed in 0(1) time assuming a data structure allowing 
0(1) time longest common prehx queries for suffixes of W and W is 
given. 


In other words, k is an efficiently-computable integer that partitions the val¬ 
ues of i into three parts: one with a single value {i = k + 1) and two others 
where either or is hxed and the other is a linear set. Handle the 

case oi i = k + 1 individually by using Lemma 4.11 to check if the word be¬ 
tween LfAkj^i and A^j^iR^^^ has a factorization into two palindromes. Next, 
check that all Ai, Ai except i = 0 are admissible by verifying L^[— 1 ] 7 ^ 

Also handle i = 0 individually, including checking admissibility. Lemma [4.9| 
is used to handle the remaining two cases in 0(log |hF|) time each. 

Case 1: 0 < i < fc. In this case, is hxed and = |Ao| -|- \R\i. 

If /c — 1 < 3, then handle all three cases individually in 0(log |hF|) time using 
Lemma 4.11 Otherwise handle only i = k similarly. 

By Lemma 4.10, A* = Aq for all 0 < i < fc — 1 and |Ao| < \R^~^\ 


Aoi?* for all 0 < i < A: — 1 and R^ are suffixes of W. Also, for all i < k 


So 

- 1 , 
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\AqR'\ < \AoR’^~^\ < \AoR^\ — |i?| < |i?^| — |i?|. So for some R' with 
\R'\ = |/?|, {Aq-R* 0 < i < k — 1} = {(-R')Mo : 0 < i < /c — 1}. Let L' be 
such that W = AiL'{R'YA q. Then a solution factorization exists if and 
only if there exist palindromes P1P2 = for some 0 < i < fc — 1. 

First, search for solutions with |Pi| > \L'\. We first prove that R! is 
primitive, allowing Lemma 4.9 to be invoked. Suppose, for the sake of con¬ 
tradiction, that R! has a period of length p with p < |P'|. So {RY has 
periods of length |P'| and p such that |P'| +p < |(P0^I so by Lemma 4.2 


has a period of length gcd(|P'|,p) < |P'|. Then since {Rf contains R as 
a subword, R also has a period of length p and thus is not primitive, a 
contradiction. 

For solutions with |Pi| < \L'\, Lemma 4.7 implies that it suffices to check 
for solutions with the longest possible Pi (longest possible P 2 is handled 
when performing the symmetric iteration over values of c). To do so, scan 
the prefix palindrome factorization starting at L'[|Pi| -|- 1] for a palindrome 
of length \L'\ — |Pi| -|- |(P')^“*| with 0 < i < A: — 1. 

Case 2: fc + 2 < i < L In this case, |AiP*| is hxed and = 

\R\ + |v4fc+2| — \R\{i — {k + 2)). Then there exists a word R such that 
IV = R(Ry-^ArRAR for all fc + 2 < i <J. Let R be the suffix of RAr 
of length \R\. Then W = RAr{L'y~^^^^^RArR^. So there exists a pair of 
palindromes Pi, P 2 with W = RAiPiP 2 AiR} for some k + 2 < i < I li and 
only if {L'Y~^R = P 1 P 2 for some k + 2 < i < 1. This situation is identical 
to that encountered in the previous case - handle in the same way. 

Handling overlap. The description of the algorithm so far has ignored 
the possibility that \RAi\ + \AiR^\ > |fF|, i.e., that RA^ and AjP* “overlap”. 
For the case of 0 <i<k, this occurs when \RAq\ + \AqR}\ > \W\. Restricting 
the values of i to satisfy 0 < i < min(A;, [(|hF| — h\L\ — 2|y4o|)/|P|J) ensures 
that W can be decomposed as claimed. For the case of A: -|- 2 < i < /, this 


occurs when IL^P” 
max(|'(2|74r 


*71^1+ |^P” 


> |fF|. Restricting the values of i to satisfy 
-I- 2r|P| -I- h\L\ — |fF|)/|P|], k + 2) < i <l ensures that W can 
be decomposed as claimed. Check the individually handled cases, namely 
i = A;, A; -|- 1, for overlap individually. 

Running time. The running time of this algorithm is 0((Z -|-r) log |hF|), 
since the amount of time spent for each value of b and c is 0(log |fF|) to 
handle individual values of i and 0(log|hF|) to handle each large case by 
Lemma 4.9 However, this assumes a data structure enabling 0(1) time 
common prefix queries on W and W. Compute such a data structure in 
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0(|W^|) time using Lemma 4.12, Since f2(|iy|) time must be spent to decide 
if a boundary word has a half-turn factorization, such a computation has no 
additional asymptotic cost. □ 


Lemma 4.14. Let P be a polyomino and B{P) = ABC AD E with B, C, D, E 
palindromes. Then the gapped mirror pair A, A and palindromes B, C, D, 
E are admissible. 


Proof. A^A is admissible. It cannot be that \B\ = \C\ = 0, since then 
ABC A = AA is non-simple. If BC\V\ = BC[—1], then \B\, \C\ > 0 and 
thus i?[—1] = ^^[l] = BC[—1] = C[l] and BC is non-simple. So BC\}] ^ 
BC[—1] and by symmetry, DE\i\ ^ DE[—1]. So A, A are admissible. 

5, C, E are admissible. Consider the pairs of non-equal consecutive 
letters in W. These pairs come from sets TZ = {lu, ur, rd, dl} and C = 
{ul. Id, dr, ru}, and Proposition 6 of [7j states that the number of pairs from 
TZ is four more than the number from C. Also, any palindrome contains an 
equal number of consecutive letter pairs from C and TZ. 

If I A| = 0, then W has factorization W = BCDE with the four consecutive- 
letter pairs from TZ not contained in any factor, i.e., for each factor X G 
{B,C,D,E], W = XV with F[-1]X[1], A:[-1]F[1] G 7Z. Since X is a 
palindrome, X[-1]X[1] = X[l]y[l] G and so X[l] 7 ^ y[-l]. Thus X is 
admissible. 

If |A| > 0, then \BC\ > 0, since otherwise A[—1]A[1] = A[—1]A[—1] is 
a subword and W is non-simple. Without loss of generality, |i?| > 0. If 
\C\ = 0, then W = BY with F[l] = A[l] = A[-l] ^ A[-l] = F[-l] and B 

is admissible. If \C\ > 0, then C[l] = C'[— 1] 7^ A[l] = A[—1]. So IT = BY 
with 1] = A[—1] 7 ^ C[l] = y[l] and so B is admissible. By symmetry, it 
is also the case that C, D and E are also admissible. □ 


Theorem 4.15. Let P be a polyomino with \B{P)\ = n. It can be decided in 
O(nlog^n) time if B{P) has a half-turn factorization. 


Proof. First, compute the prehx palindrome factorizations of each shift of 
W by computing and truncating the prehx palindrome factorizations of the 
|hF| longest suffixes of W using Lemma 4.5 Similarly compute the suffix 


By Lemma 4.5, this takes 


palindrome factorization of every shift of W. 

0 (|1F| log |hF|) total time. 

Next, compute the admissible factors (including zero-length factors), i.e., 
palindromes maximal about their center, by computing and truncating the 
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maximal palindromes of WW output by Manacher’s (9(|iy |)-time algorithm 
Each admissible factor F is contained in a prehx palindrome factorization as 
... Xl with either 1 < i < m and 0 < j < r*, or i = m and j = 
if F is the longest prehx palindrome of the word. If j < call the 
terminator of F, otherwise call the terminator of F. Either all or none 
of the prehx palindromes with a given terminator are admissible. Similar 
dehnitions and observations apply to suffix palindrome factorizations. 

For each admissible factor W, mark the two terminators (one prehx, one 
suffix) of the factor. Locating the prehx and suffix terminators for each of 
the 2|iy| admissible factors takes (9(|iy | log IhE|) total time. 

Without loss of generality, every solution half-turn factorization has |E| > 
0. Search for half-turn factorizations ABCADE by iterating over possible 
hrst letters of E. By Lemma |4.14[ only solutions with admissible D and E 
must be considered. This corresponds to a palindromes D and E starting 
and ending with a marked terminators X ^” an d Xp^, respectively. For each 
such terminator pair X_ 


■y'f'p 
) -^p 1 


use Lemma 


4.13 


with L = X, 


pi 


I Tp, R Xg 


r = to check for a partial factorization ABC A to complete the factorization 
along with D and E with the marked terminator pair. By Lemma 4.4, Xg 
and Xp are primitive and by dehnition of palindrome factorizations, I and r 
are maximal. 

Checking for a partial factorization using Lemma 4.13, each pair Xl%Xp^ 
takes 0{{rg -|-r^) log |hL|) time, 0(log|hL|) time for each admissible factor 
involved. Moreover, each admissible factor is involved in 0(log|iy|) pairs 
of terminators: 0(log|lE|) prehx (suffix) terminators when E (D). So 
0(log^ \ W\) total time is spent per admissible factor and in total the the 
algorithm takes 0(|1E| log^ |1E|) time. □ 


5 Quarter-Turn Factorizations 

Definition A quarter-turn factorization of a boundary word W has the form 
W = ABC with A a palindrome and B, C 90-dromes. 

5.1 Long 90-dromes 

Lemma 5.1. Let W be a word with a period of length p, and X a 90-drome 
subword of W. Then |X| <p. 
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Proof. Suppose, for the sake of contradiction and without loss of generality, 
that |X| = p + 2. Since W has a period of length p, X has a period of 
length p. So X[l] = X[p+1] = X[—2] and X[2] = X[2 +p] = X[—1], Since 
X is a 90-drome, X[l] = t 9 o(X[-l]) and X[ 2 ] = t 9 o(X[- 2 ]). So X[l] = 
t 9 o(X[-l]) = tgo{X[2]) = ti8o(X[-2]) = ti 8 o(X[l]), a contradiction. □ 

Lemma 5.2. Let W be a word. Let P ^pre W, S W be distinct 90- 
dromes with |P|, [S'! > 2/3|iy|. Then any other 90-drome factor of W with 
center not shared with P or S has length less than 2\W\ — (|P| -I- 151). 

Proof. Let P = PlPr and S = SlSr with \Pl\ = |Pr| and 15^1 = 15/^1. So 
there exists A such that W = PlASr (see Figure]^. Since |P/| -|- |y4| -|- 
\Sr\ = \W\, |A| < |hF|/3. Then because |Pr|, \Sl\ < |kF|/3, it follows that 
APPr.Sl. 


w 

Pl 

Pr 


A 


Sl 

Sr 


Figure 5: The words used in the proof of Lemma 5.2 


Let 1 < /,r < |hF| such that r = I 2\A\, I < \Pl\, and (/ -f r )/2 7 ^ 
(|P| — l)/2 (the center of fF[/..r] is not the center of P). The approach is to 
prove W[l] 7 ^ tgo{W[r]) and thus any 90-drome factor with center not shared 
with P, or 5 by symmetry, has length at most r — I = 21^41 = 2|fF| — (|P| -1- 
|5|) < 2|hF| — 2/3(|P| -I- |5|). Since for any pair 1 < /,r < |fF| such that 
r — I = 21^41, either I < \Pl\ or \PlA\ < r, there are two cases: I < \Pl\ — |A| 
and \Pl\ — |A| < / < |Pl|. 

Case 1: I < \Pl\ -\A\. 

W[l] = P^[l] 

= ^9o(-Pr[|-Pr| ~ 1 ]) 

= tgo{SR[\PR\-l-\A\]) 

= Sl[\Sl\-{\Pr\- 1 -\A\)] 

= P,A[\P,A\ - {\P^\ - I -\A\)] 

= PlA[1 + 2\A\] 

= hF[r] 
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Case 2: \Pl\ - |A| < / < \Pl\. 

W[l] = Pl[1] 

= ho{PR[\PR\ - /]) 

= ^9o(^[|-Pr| ~ ^]) 

= t9o(^L[|^L|-|A| + (|Pi?|-0]) 

= t90(^L[|^L| + (|Pfl|-|A|-/)]) 

= ^i8o(5'_r[—(|Pr| — \ A \ — /)]) 

= hsom\PL\ + \A\ - {\Pn\ -\A\ - 1)]) 
= tiso{W[l + 2\A\]) 


□ 

Lemma 5.3. Let P,Q,W be 90-dromes with P,Q ^pre W and \P I < IQI < 
|VL|. Then \P\ < 2/3|lL|. 


Proof. If|Q| < 2/31ly I, then clearly the lemma holds. Otherwise, Lemma 5.2 
implies that \P\ < 2|iy| - 2/3(|P| + |iy|) = 4/3|iy| - 2/3|P| < 4/3|iy| - 
2/3|iy| = 2/3|iy|. □ 


Lemma 5.4. Let W be a boundary word. There exist 0(1) admissible 90- 
drome factors of W with length at least |iy|/3. 


Proof. Consider the set of all such factors with centers contained in a 
factor of W with length at most |iy|/16. Let G be the shortest factor of 
W that contains all factors in ^. Either there exist two distinct factors 
P,S E such that P ^pre G and S ^suff G or G G J^. 

First, assume the former. Since the centers of P and S lie in a common 
factor of length at most |iy|/16, P ^pre G, and S ^sus G, |P| + |iy|/8 > 151. 
So |G| < |P|/2 + |iy|/16 + |5|/2 < |P| + |iy|/16 + (|P| + |iy|/8)/2 < |P| + 
|iy|/8 < 3/2|P| and thus 2/3|G| < |P|. By symmetry, also 2/3|G| < |5|. 
Moreover, since every element of ^ is admissible, no pair of elements share 
the same center. Then by Lemma 5^, any factor F & ^ with F ^ P,S 
has length less than 2|G| — (|P| + |5|) < |P| + |iy|/4 + |5| — (|P| + |5|) = 
|iy|/4 < |iy|/3. So ^ = {P,5}. 

Next, assume the latter. Let G' be the shortest factor with the same 
center as G' that contains all elements of except G. Since G' is not 
admissible, G' ^ J^. Without loss of generality, there exists P' G such that 
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P' ^pre G'- Using the same argument as before, 2/3|G'| < \P'\. Also, clearly 
2/31G'I < |G'| and since P' is admissible and G' is not, P' ^ G'. Then any 
factor F e with F ^ P', G', G has length less than 2|G'| — (|P'| + |G'|) = 
|G'| - \P'\ < \W\/A. So ^ = {P',G}. 

Thus the maximum number of admissible 90-drome factors of W with 
length at least |lU|/3 is at most 17|J^| = 34, obtained by partitioning all 
such factors into at most 17 groups with centers contained in common factors 
of length at most |iy|/16. □ 


5.2 Long palindromes 


Lemma 5.5. Let W be a word. Let P ^pre W, S ^sufr hh with |P|, 151 > 
2/31 IT I, P 7 ^ S, and P, S palindromes. Then W has a period of length 
2\W\-i\P\ + \S\). 


Proof. Let P = PlPr and S = SlSr with \Lp\ = \Rp\ and \Ls\ = |Ps|- The 
proof proceeds identically to that of Lemma 5^, except that since elements 
equidistant from the center of a palindrome are equal, rather than at a 90° 
rotation as in a 90-drome, W[i\ = W[i + 2|A|] for all i with either 1 < z < 
\Pl\ — |A| (case 1) or \Pi,\ — |A| < z < \Pl\ (case 2). Thus, by symmetry, 
W[i] = W[i + 2|A|] for all i with 1 < z < |1T| - 2|A|. □ 


Lemma 5.6. Let W he a word. Let P ^pre W, S ^sus hh be distinct 
palindromes with |P|, |5| > 2/3|lT|. Then any admissible middle palindrome 
factor ofW has length at most 2|1T| — (|P| -f |5|). 


Proof. By Lemma 5.5, W has a period of length p = 2|1T| — (|P| -|- |5|). Let 
1T[/ — l..r -|- 1] be a middle palindrome factor of W with length more than 
p. So (r — 1) — (/ -I- 1) -I- 1 > p and thus r — I > p + 2. 

Consider W[r—pi], W[l+pi] for integers z such that 1 < /-|-pz, r—pi < |1T|. 
Since W has a period of length p, W[l + pi] is equal for all such values of z, 
as is W[r — pi]. 

Let /(z) = (r — pi) — {I + pi). Then /(O) = r — l>p + 2 and /(z -|- 1) = 
/(z) —2p, so there exists an integer a; > 0 such that |/(a;)| < p. Also, W[l+px] 
and W[r — px] are equidistant from the center of a palindrome of length at 
least p -|- 1. So 1T[Z -|- px] = W[r — px]. Combining this with the previous 
observation, W\l] = W\l + px] = W[r — px] =W[r]. So PU[/ -|- l..r — 1] is not 
admissible. □ 
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Lemma 5.7. Let W be a word. There exists a 0{l)-sized set ^ of factors of 
W such that every admissible palindrome factor with length at least |PL|/3 is 
an affix factor of an element of ^. Moreover, ^ can be computed inO{\W\) 
time. 


Proof. Three factors. Let P, Q, S be distinct admissible palindrome factors 
of W with length at least |hL|/3 and centers contained in a factor of W with 
length at most |hL|/16. 

Let G ^ IT be the shortest factor such that P,Q,S ^ G. Now it is 
proved that if G 7 ^ P, Q, S, then P, Q, S G. 

Without loss of generality, suppose P ^pre G, S ^sufr G. Since the 
centers of P and S lie in a common factor of length at most |1T|/16 and 
S ^ G, |P| + |lT|/8 > |^|. So |G| < |P|/2 + |1T|/16 + \S\/2 < |P|/2 + 
|lT|/16 + (|P| + |lT|/8)/2 = |P| + |lT|/8 < 3/2|P| and thus 2/3|G| < |P|. By 
symmetry, also 2/3|GI < \S\. Then by Lemma 5.6[ IQI < 2|G| — (|P| + |S'|) < 
|P| + |lT|/8 + 1^1 - (|P| + 1^1) = |lT|/8. So Q G and thus Q G. 

More than three factors. Consider a set = {Pi, P 2 ,..., P^} of 
at least three admissible factors of W of length at least |lT|/3 such that 
the centers of the factors are contained in a common factor of W of length 
I IT 1/16. Now it is proved that every element of is an affix factor of one 
of the two factors. 

Let G ^ IT be the shortest factor such that PPG for every P G .T. 
Either there exist distinct P,SeJ^ with P ^pre G, S ^sus G, or that G E 
and every F E besides G has P ^mid G. 

In the first case, the previous claim regarding three factors implies P 
G for any F E JL with F ^ P,S. Also P, S ^afr G. So every factor in is 
an affix factor of G. 

In the second case, let G' ^ G be the shortest factor with the same center 
as G such that every factor in excluding G is a factor of G'. Since G' is 
not admissible, G' ^ J^. Without loss of generality, there exists P E ^ 
such that P ^pre G'. Since P E ^ and G' ^ ^, P ^ G'. Using the same 
argument as before, 2/3|G'| < |P|. Then by Lemma 5.6, every middle factor 
of G' in has length at most 2|G'| — (|G'| + |P|) < |G'| — |P| < |lT|/3. So 
every factor of G' in ^ is an affix factor of G'. Thus every factor in ^ is 
either G or an affix factor of G'. 

All factors. Partition IT into 17 factors G, G, ..., In each of length at 
most I IT 1/16. Let be the set of factors with centers containing letters 
in Jj. Then by the previous claim regarding more than three factors, there 
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exists a set [G and possibly G') such that every element of Gi is an affix 
factor of an element of and \^i\ < 2. So every admissible F ^mir fh with 
l-^l > |W^I/3 is an affix factor of an element of ^ = IJili \^\ < 2-17. 

Computing ^. The proof as described nearly yields the algorithm. Use 
Manacher’s algorithm [22] to compute all 0(|1U|) admissible palindrome fac¬ 
tors of W in 0(|1U|) time. Remove all factors shorter than |lU|/3 and group 
the remainder into 17 groups Ji,... ,/i 7 as previously described. For each 
group of palindromes, sort the indices of the hrst and last letters separately. 
Use the hrst two values in the sorted lists, including ties, to compute and 
output G, and G', if necessary. □ 


5.3 Algorithm 

Lemma 5.8. Let P be a polyomino and B{P) = ABC with A a palindrome 
and B, C 90-dromes. Then A, B, C are admissible. 

Proof. A is admissible. Without loss of generality, |i?| > 0 and either 
|G| = 0 or |G| > 0. If |G| = 0, then U[l] = R[l] = t 9 o(^[-l]) ^ to(>"[-l]). 
If |G| > 0, then U[l] = R[l] = t9o(5[-l]) ^ t 9 o(C'[l]) = t 27 o((t 9 o(C'[-l]))) = 

(o(y|-i]). 

B, C are admissible. If 1^41 = 0, then |R|, |G| > 0. So G[l] = 
^90 (^*[”1]) 7 ^ ^ 27 o(C'[—1]). So fU = BC with G[l] ^ f 27 o(C'[~l]) and thus B is 
admissible. If |R| > 0, then G[l] = t_ 9 o(G[—1]) ^ t_ 9 o(R[l]) = t 27 o(^[l]) = 
^ 27 o(^[—!])• So fU = BY with U[l] = G[l] 7 ^ t 27 o( 2 l[—1]) and thus B is 
admissible. Combining both cases, B is always admissible. By symmetry, G 
is also always admissible. □ 


Theorem 5.9. Let P be a polyomino with \B{P)\ = n. It can be decided in 
0{n) time if B{P) has a quarter-turn factorization. 


Proof. Computing admissible factors. Use Lemma 4.12 to preprocess 
W and W in 0(|hF|) time to allow 0(l)-time computation of the longest 
palindrome for each center W[i\ or W[i..i + 1]^ i.e., the admissible palindrome 
factors of W. Do the same with W and tgQ{W) to compute the admissible 
90-dromes of W. For each letter W[i\, construct a length-sorted list of the 
admissible 90-dromes that start at W[i]] do the same for those that end at 
W[i]. Also construct similar lists of the admissible palindromes. 

Long factors. By pigeonhole, any quarter-turn factorization has at least 
one long factor of length at least |hF|/3. Lemma 5.4 states that there are 0(1) 
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long 90-drome factors, while Lemma |5.3 implies that the long palindrome 
factors can be summarized by a 0(|hL|)-time-computable, 0(l)-sized set of 
letters that contains either the hrst or last letter of every such factor. Finally, 
Lemma [53] implies that for any letter, there are 0(log \ W\) 90-drome factors 
that start (or end) at the letter, and thus 0(log^ |fF|) double 90-drome factors 
that start (or end) at the letter. The approach is to search for factorizations 
containing long palindrome or 90-drome factors separately by guessing the 
90-drome factors of the factorization, given either a long 90-drome factor 
or the location of the hrst or last letter of a long palindrome factor. The 
correctness follows from observing that every factorization can be constructed 
in this way. 

Long 90-drome factors. Scan through the 0(|fF|) admissible 90-drome 
factors and extract the 0(1) long factors. For each such factor fF[L.j], use 
the previously computed lists to combine this factor with the admissible 
90-drome factors that end at W[i — 1] and start at W[i + 1], yielding all 
factorizations W = AD 1 D 2 such that Di or D 2 equals W[i..j] and both are 
admissible 90-dromes. Then check for each factorization if ^4 is a palindrome, 
again using the previously computed lists. There are 0(1) choices for W[i..j], 
0(log|hF|) factorizations containing it, and checking if A is a palindrome 
takes 0(log \W\) time. So 0(log^ \W\) total time is spent. 

Large palindrome factors. Let W[i\ be a letter in the aforementioned 
set of letters summarizing the long palindrome factors. Iterate through the 
double admissible 90-dromes D 1 D 2 starting at W[i + 1] using the previous 
lists, including choices with IO 2 I = 0. Each choice induces a factorization 
W = ADiD 2 'i check if A is a palindrome. Repeat the same process, but 
for factors D 1 D 2 that end at W[i — 1], There are 0(log^ |hF|) choices for 
the double admissible 90-domes, and checking if A is a palindrome takes 
0(log |IF|) time. So 0(log^ |IF|) total time is spent. □ 


6 Type-1 Reflection Factorizations 

Definition A type-1 reflection factorization of a boundary word W has the 
form W = ABfe{B)ACfq>{C) for some 0, <h. 
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6.1 Enumerating square and reflect square factors 

This section starts by reviewing the 0(|hh| log |lV|)-time algorithm of Main 
and Lorentz for enumerating all square factors of a word. Their algorithm 
is then modihed to output a set of 0(|1T| log \ W\) reflect square factors of 
a boundary word W that contains all factors Bf^i^B) and Cfq,{C) found in 
type-1 reflection factorizations of W. 

Let W = UV with \U\ = \V\ ± 1. The algorithm by Main and Lorentz 
hrst enumerates all squares in W that do not lie in U or V, then recursively 
calls the algorithm on U and V. The key is a 0(|hL|)-time algorithm for 
enumerating the squares that do not lie in U or V. This algorithm relies on 
the following algorithmic and structural lemmas: 

Lemma 6.1 ([25]). Let A, B be words. The longest common prefix of A and 
i?[L.|i?|] for all 1 < i < \B\ can be computed in 0{\B\) total time. 

Lemma 6.2 ([25]). Let W = UV. Let l,k be integers with 1 < / < |K| and 
I <k <21. Let Si be the length of the longest common suffix ofU and 
and Pi be the length of the longest common prefix of V and l/[i.|l/|]. Then 
the factor F with |F| = 21 and last letter at V[k] is a sguare if and only if 
21 - Si < k < I +PI+1. 


u 

V 



F[; + 1..|F|] 

V[l..k] 


Fi 

F 2 Fi 

Fi 


Figure 6: The words used in the proof of Lemma 6.2 


Proof. The factor F is a square if and only if F = FiF 2 F[F 2 such that 
Fi = F(, F 2 = F 2 , |Fi| = 21 — k, and IF 2 I = k — I (see Figure [^. In this 
case, Fi is a suffix of U, F[ is a suffix of UV[l..l], F 2 is a prehx of V, and 
F 2 is a prehx of -|- 1]. Thus the condition holds if and only if s; > 2/ — fc 
(ensuring Fi = F^') and pi+i > k — I (ensuring F 2 = F 2 ). Combining these, 
21 - Si < k < I +P1+1. □ 

Lemma 6.3 ([2S])- Let W be a word. The sguare factors of W can be 
enumerated in 0(|IF| log |IF|) time. 


25 





















Proof. A square factor F is new provided F ^ U,V . A new square factor 


of F is left provided F[l..|F|/2] ^ U and right otherwise. Then Lemma 6.2 


implies the new right square factors of length I are those ending at V[k\ for k 
in the interval [2/ — I +p/+i]. Use Lemma 6.1 to compute Pi and Sj for all 
values of i in 0(|U|) time. Enumerate the new right square factors of each 
length 21 as a set dehned by the interval [2/ — si, / + p;_|_i] of last-letter indices 
in 0(|U|) total time. 

By symmetry, computing the left square factors is possible in 0{\U\) 
time. All new square factors are either left or right, so computing all such 
factors is possible in 0{\U\ + |U|) = 0(|1U|) time. Recursion on U and V 
suffices to compute all factors in total time T(|1U|) = 2T(|lU|/2) -|-0(|1U|) = 
0(|1U| log |1U|) time. □ 


Lemma 6.4. Let W he a word. The reflect square factors of W can be 
enumerated in 0(|1U| log |1U|) time. 


Proof. For now, consider enumerating all reflect square factors XfQ{X) for 
hxed 0. Recall the proof of Lemma [6.3[ Redehne Sj in the proof to be the 
longest common suffix of fe{U) and P[l..i], and similarly pi to be the longest 
common prehx of feiP) and P[i..\P\]. Lemma 6.1 still implies pi, s* for all 
values of i can be computed in 0(|U|) total time, and the remainder of the 
proof implies all reflect square factors Xf^lX) of W can be enumerated in 
0(|1U| log |1U|) time. Perform the algorithm once for each choice of 0 to 
enumerate all reflect square factors in 0(|1U| log |1U|) time. □ 


Notice that the number of reflect square factors output is not necessarily 
0(|1U| log |1U|), and in fact may be fl(|lUp), e.g., W = (urdr)*. The running 
time is achieved by outputting factors in sets, each consisting of factors with 
the same length and contiguous last letters. Next, we prove that only factors 
of singleton set can be reflect square factors Bf^i^B) or Cfi^{C) in type-1 
reflection factorization. 


Lemma 6.5. Let P be a polyomino and B{P) = ABfQi{B)ACfq>{C). Then 
the gapped mirror pair A, A and reflect squares Bfe^B), C/$(C') are admis¬ 
sible. 


Proof. A, A is admissible. Consider i?[l] and fe{B)[—l]. Since i?[— 1] 
and fe{B)[l] are adjacent, i?[—1] feiB)[l\ and thus fe{B)[—l] .B[l]- 

Similarly, /$(C')[—1] 7^ ^[1], and thus AflA is admissible. 
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Bf(s{B), Cf^{C) are admissible. Consider the admissibility of BfQ{B) 
the other case follows by symmetry. Either 1^41 > 0 or 1^41 = 0. In the hrst 
case, suppose without loss of generality and for the sake contradiction, that 
A[l] = i?[l] = X. Then v4[—1] = 74 [ 1 ] = x and B{P) has a subword xx, a 
contradiction. 

In the second card, consider the pairs of non-equal consecutive letters in 
W . These pairs come from sets TZ = {lu, ur, rd, dl} and C = {ul. Id, dr, ru}, 
and Proposition 6 of [7] states that the number of pairs from TZ is four more 
than the number from C. Also, the total number of pairs from C and TZ in 
a pair of words X, /e(A) is equal. Thus every pair not in a factor of the 
factorization must be in TZ. 

Suppose, for the sake of contradiction, that fe{B)[l] = x and ^[l] = 
fe{x). Let B[-l] = y, so fe{B)[-l] = fe{y)- Then yx and feiyx) are 
consecutive pairs in W not in any factor (they straddle factors B, fe{B) and 
fe{B),C, respectively). Since fe{y)fe{x) is the reflection of yx, these pairs 
cannot both in TZ, a contradiction. □ 

Lemma 6 . 6 . Let W he a boundary word. A 0(|hP| log \W\)-sized superset of 
all admissible factors X = i?/©(i?) ofW can be enumerated in 0{\W\ log \W\) 
time. 


Proof. The algorithm is identical to that of Lemma 6.4, except only return 
a set if 2/ — s; = Z +pi^i, (see the proof of Lemma 6.4) i.e., there is only one 
factor in the set. All that remains is to prove that if there is more than factor 
in a set, then no factor in the set can be in a type -1 reflection factorization. 

Consider two factors lT[/..r] = Y f^fV), hP[Z-|-l..r-|-l] = Zfe{Z). Let a: = 
W[r + 1] and thus /©(F)[1] = Z[—l] = /©(x). So W[l..r] is not admissible. 


By symmetry, W[l + l..r + 1] is also not admissible. Then by Lemma 6.5 
neither factor can be in a type-1 reflection factorization. So since each set 
of reflect square factors enumerated by the algorithm of Lemma 6 ^ share a 
common length and contiguous set of last letters, the elements of the set can 
be factors of a type -1 reflection factorization only if the set is singleton. □ 


6.2 Algorithm 

Theorem 6.7. Let P be a polyomino with \B{P)\ = n. It can he decided in 
0{n\ogn) time if B{P) has a type-1 reflection factorization. 
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Proof. First, compute a set of 0(|fF| log |fF|) candidates for the admissi¬ 
ble fact ors B fejE) in 0(|fF| log |fF|) time using Lemma 6.6 Next, use 
Lemma 


4.12 


to preprocess IF, IF, and /e(lF) for 0 e {—45°, 0°, 45°, 90°} to 
enable the following queries in 0(1) time: “For Y, Z E W, W, /e(lF), what is 
longest common prehx of F[L.|F|] and Z[j..\Z\]7”. For each candidate factor 
Bfe{B), use a query on IF and IF to compute the induced A and A factors, 
and then queries on IF and /<i>(lF) for all $ to check if the remainder of IF 
is a reflection square 0/#(0). 

By Lemma 6.6, computing the candidates for i?/©(i?) takes 0(|1F| log |1F|) 
time. For each candidate, hve 0(l)-time queries are used to determine if the 
candidate can be completed into a factorization. So the algorithm takes 
0(|fF| log |1F|) total time. □ 


7 Type-2 Reflection Factorizations 

Definition A type-2 reflection factorization of a boundary word IF has the 
form IF = ABCAfQ{C)fe{B) for some 0. 

7.1 Gapped Reflect Squares 

Lemma 7.1. Let P be a polyomino and B{P) = ABCAfe{C) fe{B). Then 
gapped reflect squares B,fe{B) and C,fe{C) are admissible. 

Proof. Consider the admissibility of 5, /©(i?); the other case follows by sym¬ 
metry. Consider the pairs of non-equal consecutive letters in IF. These pairs 
come from sets TZ = {In, nr, rd, dl} and C = {ul. Id, dr, ru}, and Proposi¬ 
tion 6 of [7] states that the number of pairs from TZ is four more than the 
number from C. Also, observe that the consecutive letter pairs in factors A, 
A, B, and /©(i?) are divided evenly between £ and 7Z. Thus it cannot be 
that |A|, \C\ = 0 and three cases remain. 

Case 1: |A|, |C| > 0. Suppose, for the sake of contradiction and 

without loss of generality, that A[l] = C)!] = x. So fe{C)[l] = /©(x) and 
A[—1] = A[l] = X. Thus IF has a subword xx, a contradiction. 

Case 2: |A1 > 0, |C| = 0. Suppose, for the sake of contradiction 

and without loss of generality, that A[l] = /©(A[l]). Let x = B[—l] and 
y = A[l]. Thus /©(5)[-l] = /©(x) and A[l] = /©(?/). So xy and fe{xy) 






are consecutive letter pairs not contained in any factor of the factorization. 
Also, xy and fei^y) are not both in TZ, a contradiction. 

Case 3: |A| = 0,|C| > 0. Suppose, for the sake of contradiction and 
without loss of generality, that i?[l] = /e(C')[l]. Let x = / 0 (i?)[—1] and 
y = B[l] and thus i?[—1] = fe{x) and ^[l] = fe{y). So xy and fe{xy) are 
consecutive letter pairs not contained in any factor of the factorization and 
are not both in IZ, a contradiction. □ 

Lemma 7.2. Let W be a word. There are 0(|iy|) admissible factor pairs 
X,fe{X) ofW such that \X\ > |iy|/6. They can be enumerated in 0{\W\) 
time. 

Proof. The idea is to search for admissible factor pairs according to the 
distance d between W[a\ = X[l] and W[a + (£\ = /0(X)[1], called the match 
distance, for each d with |iy|/6 < d < 5|iy|/6. Because the factor pairs have 
|-A| > |hh|/6, any factor pair has X containing W[f\ for some t = i|iy|/6 
with i G {1, 2, ..., 6}, called an outpost index. Fixing a match distance d and 
outpost index t uniquely determines an admissible factor pair dehned by the 
largest interval W[l..r] with I <t <r such that W[j] = fe(W[j + d]) for all 
l<j<r. Thus there are 0(|fF|) such factor pairs. 

Enumerating pairs is done in 0(|fF|) time using two steps. First, prepro¬ 
cess W , /e(hF), W , and /e(hF) using Lemma 
form: 

• “What is the longest common prehx of hF[f..|hF|] and fe{W)[t+d..\W\]T'’ 

• “What is the longest common suffix of hF[l..f] and feiyV)[l..t + d]T 

Preprocessing takes 0(|hF|) total time and enables 0(l)-time queries. 

Next, iterate through all combinations of match distance d and outpost 
index t, and for each compute l,r using the data structure from the previous 
step. Output the resulting admissible factor pair X = W\l..r], fe{X) = 
W[l + d..r -I- d] unless no interval exists, i.e., W[t] ^ fe{W[t + d]). The 
0 (|fF|) combinations of d and t take 0 ( 1 ) time each and 0 (|hF|) total time 
to compute. □ 

7.2 Long admissible factors 

Lemma 7.3. Let W be a boundary word with a factor X. Let P, S W 
with such that P ^pre X, S ^sus -A, and P ^ S. Then X has a period of 
length 2|X| - (|P| + |^|). 


4.12 


to answer queries of the 
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Proof. Since P and S are mirror, there exists X' ^ W with \X'\ = |X|, 
P ^pre X', and S ^sug X'. Observe that X has a period of length r > 1 if 
and only if X[i] = X[i + r] for all 1 < i < |X| —r. Let 1 < i < |P| + [S'! — |X|. 
Then 1 < |P| + 1 - i < |X| and 1 < |P| + 1 + \S\ - |X'| - i < |^|. So: 

X[*] = P[z] 

= P[|P| + l-z] 

= JP[\P\ + l-z] 

= ^[|P| + 1 + |^|-|X'|-*] 

= ^[|^| + 1_(, + |X'|-|P|)] 

= ^[, + |X'|-|^|] 

= X[^+|X'|-|^I + (I^|-|5|)] 

= X[^ + 2|X|-(|P| + |^|)] 

Since P ^ 5, 2|X| - (|P| + |5|) > 2|X| - (2|X| - 1) = 1. So X has a period 
of length 2|X| - (|P| + |5|). □ 


Lemma 7.4. Let W be a boundary word with X P W. Let P,S P^ir fh 
such that P Ppre X, S Pguff X, and P ^ S. Any factor Y Pmid ^ with 
|F| > 2|X| — (|P| + |S'|) is not an admissible factor of W. 


Proof. By Lemma 7.3, X has a period of length r = 2|iy| — (|P| + |S'|). Let 
y ^mid X and |y| > r. 

Let X' P hh with |X'| = |X| and the center of X' exactly |hL|/2 letters 
from the center of X. Then P Ppre X', S Pgug X', and Y Pmid X'- Again 
by Lemma X' has a period of length r. 

Let U,V Y W such that W = YUYV. Since P is a middle factor of 
X, P[l..l] is in X. Since X has a period of length r and |P| > r, P[l] = 

X[|P| + 1 — r] =Y[r]. Since P is a middle factor of X' and X' has a period 
of length r, P[—1] = Y[r]. So P[l] = P[—1] and Y is not admissible. □ 


The proof of the following result is nearly identical to that of Lemma 5.7 


due to similar prerequisite structural results of Lemma 7.4 and 5.6 


Lemma 7.5. Let W be a boundary word. There exists a set ^ of 0(1) 
factors of W such that every admissible F Pmg hh with |P| > |hL|/6 is an 
affix factor of an element of ^. The set ^ can be enumerated in 0(|hL|) 
time. 
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Proof. Three factors. Let Pi, P 2 , P 3 ^mir ^ be admissible with |Pi|, IP 2 I, IP 3 I > 
I IT 1/6 and centers contained in a factor of W with length at most |1T|/14. 

Let X P IT be the shortest factor snch that Pi,P 2 ,P 3 ^ X. For some 
i,j G {1,2,3}, Pi ^pre X and Pj X. We prove that if i 7 ^ j, then 
-Pi) -^ 2 , P 3 ^aff X- 

Withont loss of generality, snppose i = 1, j = 2 and so P 3 ^mid By 
Lemma 7.4, since P 3 is admissible, IP 3 I < 2|X| — (|Pi| + IP 2 I) < |Pi| + |lT|/7+ 

IP 2 I — (|Pi| + IP 2 I) = |hb|/7 < |lT|/ 6 , a contradiction. So P 3 ^afi X. 

More than three factors. Consider a set .T = (Pi, P 2 ,..., Pm} of 
at least three admissible factors of IT of length at least |lT |/6 snch that 
the centers of the factors are contained in a common factor of IT of length 
I IT 1/14. We will prove that every element of is an affix factor of one of 
two factors of IT. 

Let G P IT be the shortest factor snch that Pj P G for every Pj G J^. 

It is either the case that there exist distinct Fi,Fr G with P; Ppre G, 

Fr Psufi G, or that G G and every Fi ^ besides G has Pj Pmid G. 

In the hrst case, Pj Pair G for any i 7 ^ /, r by the previons claim regarding 
three factors. Also P;, P^ Pair G. So every factor in ^ is an affix factor of G. 

In the second case, let G' P G be the shortest factor with the same center 
as G snch that every factor in ^ exclnding G is a factor of G'. Clearly 
G' Pmir fT and G' is not admissible. Withont loss of generality, there exists 
Fp E snch that Fp Ppre G'. Since Fp is admissible and G' is not, Fp 7 ^ G'. 

Applying Lemma 7.4 with X = G', P = Fp, S = G', every middle factor of 
G'in has length at most 2|G'| — (|G'| + |Pp|) < |G'| —|Pp| < |lT|/7 < |lT|/ 6 . 

So every factor of G' in is an affix factor of G'. Thns every factor in is 
either G or an affix factor of G'. 

All factors. Partition IT into 15 factors G, G,..., /15 each of length at 
most I IT 1/14. Let be the set of factors with centers containing letters 
in Jj. Then by the previons claim regarding more than three factors, there 
exists a set #j (G and possibly G') snch that every element of is an affix 
factor of an element of #j and |#j| < 2 . So every admissible P Pmir fT with 
l-^l P |lT |/6 is an affix factor of an element of ^ ^j and \^\ < 2-15. 

Efficient enumeration. Use Manacher’s algorithm m to compute all 
admissible (also called maximal) palindrome factors, eliminating those with 
length less than |lT|/ 6 . Partition these factors into the 15 sets G,... ,G 5 - 
For each set, sort the factors by both first and last letter, compute G and (if 
defined) G', and output them. □ 
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7.3 Algorithm 

Lemma 7.6. Let P be a polyomino and B{P) = ABCAfe{C) fe{B). Then 
A, A are admissible. 


Proof. Because of the factorization’s structure, A and A are mirror. Sup¬ 
pose, for the sake of contradiction and without loss of generality, that B[l] = 
CFn_Let ^ = femi]. So fe{C)[-l] = /e(C[-l]) = feiW]) = 
= X. Thus xx is a subword of B{P), a contradiction. □ 


Theorem 7.7. Let P be a polyomino with \B{P)\ = n. It can be decided in 
0(n) time if B{P) has a type-2 reflection factorization. 


7.6 


Proof. Recall that a type-2 reflection factorization of W has the form W = 
ABCAfe{C)fs{B). Without loss of generality, either \A\ > |iy|/6 or \B\ > 
\W\/6 . 

Case 1: |A| > |lT|/6. The factors A and A are admissible by Lemma 
Observe that each admissible factor has a distinct center iy[i..j] with j G 
{i, i-fl}. Moreover, it is the longest factor X = LR with this center s uch th at 
the factor X = RL has center W[i + \W\/2..i -|- |lT|/2]. Use Lemma 4.12 to 
preprocess VL, fU, fU, and W in 0(|fU|) time to then compute the admissible 
factor with each center in 0(|1U|) total time. 

Next, use Lemma to compute ^, a set of factors such that every factor 
A with 1^41 > |lU|/6 is an affix factor of element of ^. This takes (9(|fU|) 
time. For each element F G ^, the prehx admissible factors of F are handled 
together in 0(|hF|) time. The suffix factors of each F are handled similarly 
and symmetrically. Since \J^\ = 0(1), 0(|kF|) total time is spent. The 
remainder of this case describes how the prehx admissible factors of a single 
factor F are handled. The key is to enumerate and process each candidate 
factor pair B, /©(R), rather than candidate factor pairs A, A. 

Let W = XY such that |U| = ±|X| and F,X start at the same letter 
of W (see Figure [^. Reuse the data structure previously constructed by 
Lemma 4.12 to support the following queries in 0(1) time: “What is the 
length of the longest common suffix of feiY) and X[l..i]?”. For each i with 
1 < i < |W|, denote the result of this query Sj. 

Each pair of longest common suffixes X[i + 1 — Si..i], ^[lUl -|-1 —Si..|U|] is a 


candidate factor pair B, /©(R). Since R, f@{B) are admissible by Lemma 7.1 
no non-longest suffix can be candidates. For each pair, verify that A = 
X[l..i — Sj], the prehx of X proceeding R, is an admissible prehx factor of F. 
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Figure 7: The words used in the proof of Theorem 7.7 


This takes 0(1) time after precomputing a look-up table of the admissible 
prehx factors of F. Since A is a mirror, the prehx of Y with length 1^41 is A. 
Finally, verify that X[i + l..|X|] = / 0 (F[|A| -|- l..|T| — Sj]), i.e., that the re¬ 
maining words before and after A are factors C and fe{C), respectively. This 


is done in 0(1) time by using the data structure constructed by Lemma 4.12 


If both verihcations succeed, then the named factors form a type-2 re¬ 
flection factorization of W. In total, 0(1) time is spent per choice of i, the 
index of W where B ends, and thus 0(|kF|) time total. 

Case 2: |.B1 > |VF|/6. Compute a 0(|kF|)-sized superset of all factor 
pairs B, fe{B) using Lemma 7.2 For each pair, repeat the verihcations done 
in case 1. Using the same analysis as case 1, this takes 0(|1U|) total time. □ 


8 Type-1 Half-Turn-Reflect ion Factorizations 

Definition Let IF be a boundary word. A type-1 half-turn-reflection fac¬ 
torization of a boundary word W has the form W = ABCADfQ[D) with B, 
C palindromes. 

Lemma 8.1. Let P be a polyomino and B{P) = ABC AD f^^D) with B, 
C palindromes. Then reflect square 0 / 0 ( 0 ), and palindromes B, C are 
admissible. 

Proof. Df&{D) is admissible. Suppose, for the sake of contradiction and 
without loss of generality, that / 0 (O)[ 1 ] = / 0 (A[ 1 ]). Let x = A[l]. Then 
A[-l] = X and 0[1] = / 0 ((/ 0 (O)[l])) = / 0 ((/ 0 (A[l]))) = A[l] = x. So xx 
is a subword of B{P), a contradiction. 
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B, C are admissible. Consider the admissibility of B; the other case 
follows by symmetry. Suppose, for the sake of contradiction, that ^[—1] = 
C[l]. Let X = ^[—1]. Then y4[l] = x and C'[—1] = C[l] = x. So xx is a 
subword of B{P), a contradiction. □ 

Theorem 8.2. Let P be a polyomino with \B{P)\ = n. It ean he decided in 
O(nlogn) time if B{P) has a type-1 half-turn-reflection factorization. 


Proof. The algorithm has two phases. Phase 1 is computing a O(|IT| log \ W\)- 
sized set of factors such that if any of these is can be factored into two ad¬ 
missible palindromes, then W has a type-1 half-turn-reflection factorization. 
Phase 2 is deciding if any of these factors has such a factorization. 

Phase 1. A type-1 half-turn-reflection factorization of W has the form 
W = ABCADf 0 {D) with B, C palindromes. The factor Df^i^D) is admis¬ 
sible by Lemma 8.1 [ Use Lemma 6.6 to enumerate a 0(|1U| log |lU|)-sized 


superset of all admissible factors DfejD) 


Preprocess W and W using Lemma 4.12 to answer the following queries in 


0(1) time: “What is the longest admissible pair A, A with A ending at W[i] 
and A starting at W[i]T. For each candidate Dfe^D) = W[l..r], compute 
the longest admissible pair A, A with A ending at W[l — 1] and A starting at 
lU[r-|-l]. This partial factorization ADfQ{D)A can be completed into a type- 
1 half-turn-reflection factorization if and only if hF[r-|-l-|- |A|../ —1 — |A|] has a 
factorization into two palindromes. Call the set of such factors computed for 
all candidate Df^i^D) completion factors. These factors take 0(|lU| log |1T|) 
total time to compute. 

Phase 2. The factors B, C of a type-1 half-turn-reflection factorization 
are admissible by Lemma 8.1 Use Manacher’s algorithm [26] to compute the 


0(|lU|)-sized set of all admissible factors. For each letter in W, create length- 
sorted lists of the admissible palindrome factors that start and end at this 
letter. Constructing these lists takes 0(|1U|) time using counting sort, since 
there are 0(|1U|) factors of integer lengths between 1 and |1U|. Construct 
similar lists for the completion factors from phase 1 in 0(|1U| log |1U|) time. 

By Lemma 4T, a completion factor has a factorization into two admissible 
palindrome factors if and only if such a factorization exists utilizing the 
longest admissible prehx or suffix palindrome of the completion factor. For 
each letter in W, use a two-hnger scan to compute the longest admissible 
prehx palindrome for each completion factor starting at this letter. For each 
pair, verify that the remainder of the factor is also a admissible palindrome 
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factor using a precomputed look-up table. If so, then a type-1 half-turn- 
reflection factorization has been found. Do a symmetric scan through the 
completion factors and admissible palindrome factors ending at the letter. 

Scans take time linear in the number of admissible palindromes and 
completion factors involved. Since each factor appears in two scans (one 
each for the letter the factor starts and ends at) the total time taken is 
0{\W\\og\W\). □ 

9 Type-2 Half-Turn-Reflection Factorizations 

Definition A type-2 half-turn-reflection factorization of a bonndary word W 
has the form W = ABCDfs{B) fq>{D) with A, C palindromes and 0° —<l>° = 
±90°. 

Lemma 9.1. Let U, V be words. Let the prefix and suffix palindrome factor¬ 
izations of U and V be given, along with a data structure supporting 0(1)- 
time longest common prefix gueries for subwords of V and feifJ). It can be 
decided in 0(log \U\ ±log |D|) time if there exist factorizations U = fs{D)A, 
V = CD with A, C palindromes. 

Proof. Solve by considering solntions snch that C ends with each of the 
0(log |D|) repeated factors Xfl in the prefix palindrome factorization of V. 
If Tj < 2, then check if either choice of C ending at one of the factors leaves 
a factor D ±sufr V snch that /$(0) is a prehx of U and the suffix of U 
with length \U\ — \D\ is a palindrome. Both checks together take 0(log \ U\) 
time by using a single longest common prehx qnery and a scan of the snffix 
palindrome factorization of Lf. 

Otherwise ri > 3. Compnte the largest k snch that the longest common 
prehx of Xfl and f^{U) is at least k\Xi\ using a longest common prehx qnery. 
If /c = 0, then no factorization V = CD exists nsing a factor D starting at a 
repetition of Xp, so assnme that k > 1. Let Q be the snffix of V following Xfl. 
By the maximality of r* that follows from the dehnition of a prehx palindrome 
factorization, either: (1) IQI > \Xi\ and Xj 2<pre Q or (2) |(5| < |Xj|. Also, in 
both cases, no factorization V = CD exists such that D has a prehx X^^^ or, 
equivalently, such that D starts at the jth repetition of Xj, for any j < ri — k. 

Case 1: |Q| > |Xj| and Xj 2^pre Q- In this case, D also cannot start 
at the jth repetition of Xj for any j > r^ — k 1 and so D must start at the 
repetition r* — /c ± 1 of Xfl. Thus the factors C, D, and f^{D) are hxed. 
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Search for a suffix factor of U with length \U\ — \D\ hj scanning the suffix 
palindrome factorization of U. 

Case 2: |Q| < Check if the longest common prehx of Q and 

f<s,{U) is Q, i.e., if /$((5) is a prehx of U. If not, then D cannot start the 
jth repetition for any j > Vi — k + 1 since /$(Xf) is a prehx of U. Thus 
D must start at repetition rj — /c + 1; search for a factor A to complete the 
factorizations in 0(log \U\) time as before. 

Otherwise Q is a prehx of Xi and thus D can start at the jth repetition 
of Xi for all j > Ti — k + 1. Check for the existence of a choice of C, D such 
that U = fip{D)A and A a palindrome by scanning the suffix palindrome 
factorization of C for A with IAI e {\U\ — {{ri—j+l)\Xi\ + \Q\) \ j > r*—A;+l}. 
Since the lengths of the palindromes for each repeated factor also form a 
linear set, this takes 0(1) per repeated factor and 0(log \ U\) total time. 

Running time. Each repeated factor except the last has either < 2 or 
case 1 applies. Moreover, each can be handled in 0(1) time plus a search in 
the suffix palindrome factorization of U for a palindrome of length |f/| — |0|. 
As i increases, the length of the suffix palindrome searched for decreases, so 
these can be handled together in 0(log|0|) total time. Thus all repeated 
factors except (possibly) the last can be handled in time 0(log \ V\ +log \ U\), 
while the optional handling of case 2 takes 0(log \ U\) time. 

□ 

Lemma 9.2. Let B, U he words. Let the prefix and suffix palindrome fae- 
torizations of U be given, along with a data structure supporting 0{l)-time 
longest common prefix gueries for B and subwords of U. It can he decided 
in 0(log \ U\) time if there exists a factorization U = ABC with A, C palin¬ 
dromes. 

Proof. Solve by considering solutions such that B starts with each of the 
0(log \ U\) repeated factors Xfi in the prehx palindrome factorization of V. 
If Tj < 2, then check if either choice of B starting at one of the repetitions 
of Xi leaves a remainder of U, called C, that is a palindrome. Both checks 
together take 0(log \ U\) time by using a single longest common prehx query 
and a scan of the suffix palindrome factorization of Lf. 

Otherwise ri > 3. Compute lcp(X[%i?), the longest common prehx of 
Xfi and B. Let k be the largest integer such that lcp(X[% B) > /i;|Xi|. Either 
\cffiXl\B) < \B\ is less than or equal to \B\. 

Case 1: \cp{Xf%B) < |R|. No solution B can start before repetition 
Tj — /c + 1, since B is not a prehx of the remaining word. Also, no solution 
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B can start after repetition r* — /c + 1 by the maximality of r*. So check 
for a solution with B starting at repetition r* — /c + 1 in 0(log \ U\) time by 
searching in 0(log \ U\) time for a suffix palindrome C of length \U\ — \AB\ 
in the suffix palindrome factorization of U. 

Case 2: \cp{X^\B) = \B\. Observe 5 is a prehx of the suffixes of U 
starting at repetitions 1, 2,..., r* — /c + 1. Moreover, B is not a prehx of any 
suffix of U starting at a repetition after r* — A; + 1 by the maximality of rj. 
The lengths of the suffixes following those choices of B starting at repetitions 
1, 2,... Tj — /c + 1 form a linear set. Check for a solution C with length in 
this linear set. This check is possible in 0(log \U\) time, since each repeated 
factor in the suffix palindrome factorization can be checked in 0(1) time by 
solving a system of two linear equations. 

Amortized 0(l)-time search for C. In both of the above cases, the 
length of C completing a solution factorization is equal to \U\ — \AB\, where 
\B\ is hxed and |A| is a prehx ending at a repetition of Xj. So the lengths of C 
strictly decrease as i increases. Maintain a pointer into the suffix palindrome 
factorization of U for searching for C across repeated factors of the prehx 
palindrome factorization of U. This pointer visits each factor once, and thus 
the searches take 0(log \U\) total time and 0(1) amortized time each. So in 
total 0(1) amortized time is spent per i, and thus 0(log \U\) time total. 

□ 

Lemma 9.3. Let P be a polyomino and B{P) = ABCDfe{B)fq>{D) with 
A, C palindromes. Then gapped reflect squares B,fe{B) and 0,/$(0), and 
palindromes A, C are admissible. 

Proof. A, C are admissible. Since A and C are symmetric, proving the 
claim for A suffices. Suppose, for the sake of contradiction and without loss 
of generality, that /$(0)[—1] = i?[l]. Let x = fe{B)[l]. Then 0[—1] = 
/#(/e(2^)) = d:. So xx is a subword of W, a contradiction. 

B,f@{B) are admissible. Suppose, for the sake of contradiction and 
without loss of generality, that /©(O)!]) = /$(0)[1]. Let y = 0[1]. Then 
U{D)[1] = fe{y) and thus D[l] = U{fe{y)) = V- Also, y = C[l] = C[-l]. 
So yy is a subword of hh, a contradiction. 

D,f^{D) are admissible. Suppose, for the sake of contradiction and 
without loss of generality, that A[l] = /$(/©(!?)[1]). So A[l] = i?[l]. Then 
if z = A[l] = A[—1], 2 :^ is a subword of W, a contradiction. □ 

Theorem 9.4. Let P be a polyomino with \B{P)\ = n. It can be decided in 
0{n\ogn) time if B{P) has a type-2 half-turn-reflection factorization. 
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Proof. Recall that a type-2 half-turn-reflection factorization of W has the 
form W = ABCDfe{B) fis,{D) with A, C palindromes and 0° — <h° = ±90°. 
Without loss of generality, an element from {A, B,C, D} has length at least 
|hh|/6. The case of \C\ > |W|/6 is symmetric to that of l^l > |hh|/6; the 
remaining three cases are each handled in 0(|1R| log \ W\) time. 

All three cases use 0(l)-time queries of the following form: “For Y, Z E 
W, /©(VF), /$(hF), W, what is the longest common prefix starting at Y[i] and 


Z[j]?”. Such queries are possible after 0(|hF|) preprocessing by Lemma 4.12 


The three cases also use prefix and suffix palindrome factorizations of sub¬ 


words of W. These can be precomputed in 0(|hF| log |hF|) time by Lemma 4.5 


The case of |A| > |lF|/6 is reducible to the subroutines used to solve the other 
two cases. For this reason, the other cases are considered first. 

\B\ > \W\/6 . The pair of factors B, fe{B) are admissible by Lemma[973 


Enumerate a 0(|fF|)-sized superset of all pairs B, /^{B) in 0(|hF|) time us¬ 
ing Lemma 7.2 Determining if a given pair B,fQ{B) can be completed into 
a type-2 half-turn-reflection factorization is equivalent to solving an instance 
of the following problem: “Given words [/, V and their prefix and suffix palin¬ 
drome factorizations, do there exist factorizations U = /$(D)A, V = CD 
with A, C palindromes?”. Each such instance is solved in 0(log \U\ ±log |F|) 
time by Lemma 9.1 Thus the entire case is handled in 0(|hF|(log |17| ± 


log |F|)) = 0(|hF| log |hF|) time. 

\D\ > |VF|/6. As in the case of \B\ > |fF|/6, enumerate a 0(|fF|)- 
sized superset of all D,f^{D) in 0(|fF|) time. Handle each pair D,fip{D) 
individually as an instance of the following problem: “Given words U, B 
and the prefix and suffix palindrome factorizations of U, does there exist a 
factorization U = ABC with A, C palindromes?”. Each such instance is 
solved in 0(log \U\) time by Lemma 9.2 Thus the entire case is handled in 
0{\W\ log \U\) = 0{\W\ log |1F|) time. 


|A| > |VF|/6. First, Use Lemma 5.7 to compute a 0(l)-sized set of 
factors of W such that every factor A is an affix factor of element of 

Next, process the prefix palindromes of each F G The first letter 
W[f] of each F defines a last letter W[f — 1] for /$(D). Gombine this last 
letter of f<s>{D) with the 0(|1U|) choices for the last letter of D to obtain 
0(|hF|) pairs D,/$(D). Gompute each pair D,f^{D) in 0(1) time using 
the previously computed suffix tree data structure. Now handle each pair in 
0(log|H7|) time exactly as done for the case of \D\ > |H7|/6. The 0(|hF|) 
pairs are handled in 0(|hF| log |hF|) time total. 
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Finally, process the suffix palindromes of each F G ^. Each -F[—1] = 
W[f] dehnes the hrst letter W[f + 1] oi B. Combine this hrst letter of B 
with the 0(|hF|) choices for the hrst letter of fe{B) to obtain 0(|hF|) pairs 
B,f 0 {B). Compute each pair B^f^^B) in 0(1) time and handle them in 
0(log |fF|) total time as done in the case of \B\ > |hF|/6. □ 

10 Conclusion 

This work demonstrates that not just polynomial, but quasilinear-time al¬ 
gorithms exist for deciding tiling properties of a polyomino. It remains to 
be seen if a linear-time algorithm exists, or whether a super-linear lower 
bound for one of the factorization forms exists. The slowest algorithm is for 
half-turn factorizations, so it seems natural to attack this special case hrst. 

Open Problem 10.1. Can it be decided in o{n\o^ n)-time if a polyomino 
P with \B{P)\ = n has a half-turn factorization? 

Open Problem 10.2. Can it be decided in 0{n)-time if a polyomino P with 
\B{P)\ = n has an isohedral tiling of the plane? 

For monohedral tilings containing only translations of the prototile, a 
polyomino has such a tiling only if it has one that is also isohedral [21 ES] ■ 
Does this remain true for tilings using other sets of transformations of the 
prototile? Modifying the anisohedral tile of Heesch [18] (see [16]) proves 
that the answer is “no” for tilings with rehected tiles, while an example of 
Rhoads [22] proves that the answer is “no” for tilings with 90° rotations of 
tiles. This leaves one possibility open: 

Open Problem 10.3. Does there exist a polyomino P that has a tiling 
containing only translations and 180° rotations of P and every such tiling is 
anisohedral? 

As mentioned in Section there are isohedral tiling types (characterized 
by boundary factorizations) that cannot be realized by polyominoes due to 
angle restrictions. Moreover, the boundary factorization forms here also ap¬ 
ply to general polygons, under appropriate dehnitions of “boundary word”. 
Extending the algorithms presented here to polygons, along with develop¬ 
ing algorithms for the remaining boundary factorizations is a natural goal. 
However, signihcant challenge remains in efficiently converting a polygon’s 
boundary into a word that can be treated with the approach used here. 


39 


Open Problem 10.4. Can it he decided in 0(nlog^n) time if a polygon 
with n vertices has an isohedral tiling of the plane? 

Observe that pairs of tiles in a tiling that can be mapped to each other via 
a symmetry of the tiling indnces a partition of the tiles. Dehne a tiling to be 
fc-isohedral if the partition has k parts, e.g., an isohedral tiling is 1-isohedral. 
Thus fc-isohedral tilings are a natural generalization of isohedral tilings that 
allow increasing complexity; specihcally, they cannot be characterized by 
a single boundary factorization. A natural generalization of the problem 
considered here is as follows: 

Open Problem 10.5. Can it be decided efficiently if a polyomino has a 
k-isohedral tiling? 

An approach described by Joseph Myers [29] achieves a running time of 
approximately \ though a precise analysis of the running time has not 
been performed. A hxed-parameter tractable algorithm also may be possible. 
On the other hand, a proof of NP-hardness is unlikely, since it implies, for 
each c G N, the existence of prototiles whose only tilings are fc-isohedral for 
k > c. Such tiles are only known to exist for c < 10 m- 
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